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ABSTRACT 

At low Eddington ratios (m), two effects make it more difficult to detect certain AGN given 
a particular set of selection methods. First, even allowing for fixed accretion physics, at low 
m AGN become less luminous relative to their hosts, diluting their emission; the magnitude 
of the dilution depends on host properties and, therefore, on luminosity and redshift. Second, 
low-f7i systems are expected and observed to transition to a radiatively inefficient state, which 
changes the SED shape and dramatically decreases the luminosity at optical through infrared 
wavelengths. The effects of dilution are unavoidable, while the precise changes in accretion 
physics at low m are somewhat uncertain, but potentially very important for our understanding 
of AGN. These effects will have different implications for samples with different selection 
criteria, and generically lead to differences in the AGN populations recovered in observed 
samples, even at fixed bolometric luminosity and after correction for obscuration. Although 
the true Eddington ratio distribution may depend strongly on mass/luminosity, this will be 
seen only in surveys robust to dilution and radiative inefficiency, such as X-ray or narrow-line 
samples; by contrast, selection effects imply that AGN in optical samples will have uniformly 
high Eddington ratios, with little dependence on luminosity, even at low Lboi where the median 
"true" m < 0.01. These same selection effects also imply that different selection criteria pick 
out systems with different hosts: as a result, the clustering of low-luminosity optical/infrared 
sources will be weaker than that of X-ray sources, and optical/IR Seyferts will reside in more 
disk-dominated galaxies, while X-ray selected Seyferts will be preferentially in early-type 
systems. Taken together, these effects can naturally explain long-standing, apparently contra- 
dictory claims in the literature regarding AGN Eddington ratio distributions, host populations, 
and clustering. Finally, we show that if current observed Eddington ratio distributions are cor- 
rect, a large fraction of low-luminosity AGN currently classified as "obscured" are in fact 
radiatively diluted and/or radiatively inefficient, not obscured by gas or dust. This is equally 
true if X-ray hardness is used as a proxy for obscuration, since radiatively inefficient SEDs 
near m ~ 0.01 are characteristically X-ray hard. These effects can explain most of the claimed 
luminosity/redshift dependence in the "obscured" AGN population, with the true obscured 
fraction as low as ~ 20%. 
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1 INTRODUCTION 

At the highest luminosities, it is generally accepted that active 
galactic nuclei (AGN) must lie within a relatively narrow range of 
black hole masses and Eddington ratios - after all, a quasar with 
Lboi ~ 10 47 ergs" 1 requires a > 10 9 Mq black hole (BH); and since 
this is already on the exponential tail of the black hole mass func- 
tion, there are vanishingly few ~ 10 Mq BHs that could give a 
similar luminosity at low (~ 0.01) Eddington ratios. At low lumi- 
nosity, however, there is an obvious degeneracy between BH mass 
and Eddington ratio - a Seyfert (Lboi < 3 x 10 45 ergs~') could be 
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a low-mass BH at high Eddington ratio, or a high-mass BH at low 
Eddington ratio. Moreover, although there are more low-mass BHs 
by number density, the duty cycle of low-E ddington ratio activity 
is higher in many physica l models (see e. g. Di Matteo et alj 20031 . 
| 2008|- iGranato et alj|2004 iHopkins et afll2005 d ; lLapi et alJhoOfj 
Isiiacki et alj|2007l) , so the two may contribute comparably or in 
different proportions as a function of luminosity or redshift. 

Although it is often implicitly assumed that Seyfert samples 
selected in the same bolometric (obscuration-corrected) luminosity 
range are reasonably homogeneous, there are at least two effects - 
one observational and one physical - that could introduce dramatic 
differences depending on the wavelengths of identification and the 
methods used to identify AGN. 
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Observationally, even if accretion were self-similar, i.e. there 
were no change in SED with m, many samples implicitly require 
that the ratio of the AGN optical continuum to host luminosity be 
above a certain threshold for a galaxy to be identified as an AGN. 
This criterion varies between different samples and wavelengths, 
both because of physical differences in the SED shapes of AGN and 
galaxies, and because of various survey depths and selection meth- 
ods. Most limited are surveys which select on the basis of "typical" 
quasar colors dominated by the optical/UV continuum: this implic- 
itly requires that the AGN continuum in the bands of interest be at 
least comparable to, if not greater than, that of the galaxy. 

But even a complete spectroscopic survey looking for Type 1 
systems (i.e. the identification of broad lines above some host con- 
tinuum) implicitly requires an AGN continuum luminosity greater 
than ~ 0. 1 times that of the host (or else the lower-level broad re- 
gions of the lines can disappear into the host continuum; see e.g . 
IVanden Berk et alj|2006l ; ISanchez et al]|2004 Ijahnke et al.ll2007l) . 
This implicitly sets an Eddington ratio cut as a function of host 
galaxy properties, where the Eddington ratio m (the dimensionless 
mass accretion rate) is defined as 
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(where the last equality comes from assuming a standard radiative 
efficiency of e r ~ 0. 1 

To assess the eff ects of dilution quantita tively, assume that 
a system lies on the iMagorrian et alj dl998l) relation between 
BH mass and host bulge stell ar mass A/bh ~ (iM..^ (where 
(i ~ 0.001; Haring & Rix 2004). The total bulge mass is related 
to the total stellar mass by the bulge-to-total stellar mass ra- 
tio B/T = A/*,bui/A/», and the total host galaxy luminosity will 
be correspondingly given by the stellar mass-to-light ratio in the 
observed band, (M»/L)h os t- The total host luminosity in the ob- 
served band is therefore Mbh/[m (M*/L)h 0s t x B/T]. Assuming for 
simplicity that the radiative efficiency is ~ 0.1, so that A = m, 
the AGN bolometric luminosity is simply L = mL^i = m(3.3 x 
10 4 Lq) (A/bh/A/q). If the AGN bolometric correction to the ob- 
served band is C\ = Lband/iboi, then the ratio of the AGN lu- 
minosity in the observed band, to the host galaxy luminosity in 
the same band, will be proportional to the Eddington ratio and 
given by ^ « 33mC x B/T (M*/L) kost (with (M„/L) hos , in so- 
lar units). For the B-band, with typical parameters C\ ~ 1/12 
and (M t /L)host ~ 0-5 — 3Mp) /Lp) (for bol ometric corrections, see 
iRichards et alj2006l :l Hopkins et al . 2007b), this translates to 
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Given typical mass-to-light ratios the AGN will only be visible in 
color selected samples or identifiable as a Type 1 system by most 
survey criteria for high accretion rates m > 0.01 — 0. 1 . Note also the 
scaling with the host mass-to-light ratio - a younger or more star- 
forming galaxy is brighter, masking more of the AGN continuum 



emission, yielding a potential bias a gainst such systems in ce rtain 
samples of optical/IR-selected AGN jjahnke et al J2004ll2006h . By 
contrast, in the X-rays galaxies have much higher M/L - a sys- 
tem with an X-ray luminosity ~ 10 43 ergs -1 will be well above the 
luminosity of even a strongly star-forming host. 

Physically, at low Eddington ratio, models predict a tran- 
sition from a standard, optically thick, UV-bright thin disk 
to a radiatively ineffic i ent, optically thin accretion flow (see 



Naravan & Yi 1994; Esin et alj[l997|; lYuan & Naravanll2004l; 



lOuataert & Narayar] 1 19991, and references ther ein), as seen in 
X-ray binaries ( McClintock & Remillard 2006); many observa- 
tions of both AGN and X-ray binaries support this scenario 
dNaravan et alJl99d;feej200ll;lHol2002l;lMaccarone et alj2003l ; 



iMarchesini et alj2004IJestej|2005l ; ICao & Xull2007ri . 

In radiatively inefficient accretion flows (RIAFs), the thin disk 
is assumed to evaporate and/or move to much larger radii, and the 
broad line region and most optical/UV emission declines, while a 
hard X-ray spectrum from the thick disk remains. In such a state, 
the system could still appear in X-ray samples with nearly the 
same hard X-ray luminosity as it would have if the radiatively ef- 
ficient solution were extrapolated to the same physical accretion 
rate, but the system would disappear from optical samples look- 
ing for optical/UV light above the host galaxy continuum or for 
broad-line emission. The systems would likely, however, still ap- 
pear in narrow-line AGN samples, since the central hard X-ray 
source remains a significant source of photoionization (exactly how 
the narrow-line properties scale for an incident RIAF spectrum is 
not certain; a detailed investiga tion of these sailings will be the 
subject of future work). Indeed, iKewlev et"al, argue that the 

distinction between Seyferts and LINERs largely owes to Edding- 
ton ratio, with LINERs being lower Eddington ratio systems with a 
harder photoionizing spectrum, consistent with our discussion here. 

The models and observations (see references above) indicate 
that such a transition occurs around an accretion rate m ~ 0.01. 
Although the role of such radiatively inefficient states remains con- 
troversial, if such a change does occur, it would clearly rule out 
such a mode contributing to the X-ray equivalents of bright quasars. 
As above, it would require enormous BHs to produce such a lumi- 
nosity at such low Eddington ratio - but RIAFs could easily con- 
tribute significantly to lower-luminosity samples. A ~ 10 s Mq BH 
is still radiating at nearly 10 43 ergs~' in hard X-rays at m = 0.01, 
and the total number density of such systems is ~ 10~ 2 Mpc -3 . 
Only ~ 1 — 5% of them need to be active at these low Eddington 
ratios to account for the entir e observed popula tion of such X-ray 
sources in the local universe dUeda et alj[2003b . But with rapidly 
declining optical luminosity and disappearing broad line regions, 
such sources may vanish from most Type 1 optical samples. 

Together, these effects could lead to very different Seyfert 
populations being observed in different bands or with different se- 
lection criteria. They also relate to a physical distinction in how 
AGN may be fueled. A low-mass BH at high Eddington ratio is 
likely to be in a gas-rich, disk dominated system (given the typical 



We also define the dimensionless Eddington-scaled luminosity, A, in 
terms of the bolometric luminosity 

\ = Z-b0[ = „. 

~ L Ed d ~ 1.3 X 10 46 ergs- 1 (M BH /1O 8 M )' 

This is equivalent to m when the radiative efficiency is constant, i.e. for m 
above the critical value ~ 0.01, A = in. At lower m, A may drop much more 
rapidly than m, if sources become radiatively inefficient (and A inferred 
from different bands will be different). 



2 Note that, in realistic models, the broad line region does not simply van- 
ish instantly, but rather moves outwards with the thin disk and so declines 
with the optical/IR luminosity of the AGN. Of course, this decline, relative 
to the host, will still yield a system that will not appear in traditional broad- 
line samples. Observationally, the presence of some weak broad lines, even 
in very low accretion rate systems (many of which show continuum or po- 
larization evidence of being in a RIAF/thick disk state inwards of the thin 
disk) has been noted down to accretion rates in ~ 0.001 (see e.g. lHo|[l999r . 
iMoran et alj200d ; lTrarl2003r . lHao et alj2005bll3) . 
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properties of the expected low bulge-mass host). Fueling growth 
at these luminosities and masses need not be related to violent 
events such as major mergers of galaxies, but may be achievable by 
commo n mechanisms such as the stochastic accretion of molecular 
clouds l Hopkins & He rnquist 2006), gas inflow in barr ed systems 
dJogeell2006ir and minor mergers JVounger et alj|200^) . Although 
the fueling mechanisms can be relatively mundane, the small-scale 
physics present are - being systems at high Eddington ratio - sim- 
ply a scaled-down version of quasar activity, and as such repre- 
sent a continuation of a feedback-regulated scenario, where high- 
Eddington ratio activity (implying rapid BH growth) is truncated by 
AGN feedback at a critical mass/luminosit y corresponding to th e 
observed MBH-host spheroid correlations (Di Matteo et alj|2005b . 
and the AGN rapidly decays to lower lumino sities with some char- 
acteristic po wer law-like light curve decline ( Hopkins & Hernquist 
l2006ll2009l) . 

On the other hand, a high-mass BH at a low Eddington ra- 
tio, living in a massive host, is a natural consequence of a sys- 
tem already in its late-time decay phase, having previously gone 
through a state of rapid accretion, feedback, and self-regulated 
growth. It is possible, of course, to re-ignite such a system by the 
accretion of new cold gas in a minor merger or via secular in- 
stabilities, but this will more likely lead to a brief new burst of 
bright quasar activity before declining again to a long-lived low 
Eddington ratio mode. Simulations and observations suggest that 
syste ms can linger in this "decaying" mode for ~Gyr to a Hubble 
time dYu et alj|2005fc IVolonteri et al|2006l : IHopkins et al.ll2005bl R 
2008b:; lYu & L j 120081 : ICiotti et aljfcofm either in such a post- 
merger steadily declining level of activity or eventually in an equi- 
librium, accreting e.g. hot gas from the diffuse, pressure supported 
gaseous halo in the spheroid or from stellar winds (the distinction 
between this "steady state" accretion and a "smooth decline" from 
initial high levels of activity triggered in e.g. mergers is largely se- 
mantic - the two likely blend seamlessly into one another). 

Given the potential importance of these observational selec- 
tion effects, and the theoretical interest in understanding differ- 
ent modes of accretion and their implications for observed AGN 
samples, it is of considerable interest to develop predictions and 
observational discriminants between these modes of accretion; it 
is also important to recognize the observational selection effects 
that might lead to the preferential selection of different popula- 
tions. Any observations intending to study AGN host galaxies, fu- 
eling mechanisms, and evolution must account for such distinc- 
tions. In this paper, we adopt a simple and general model for AGN 
lightcurves and use it to illustrate the effects of dilution and radia- 
tive efficiency on observational inferences of the Eddington ratio 
distributions, obscured fractions, clustering, and host galaxy mor- 
phologies of AGN as a function of luminosity and redshift. Recog- 
nizing these biases, however, the same observations, coupled with 
models such as those presented here, can be used to break degen- 
eracies between e.g. populations with different BH masses and Ed- 
dington ratios at fixed luminosity; these can ultimately be used to 
constrain models for the long-term evolution of AGN across many 
different stages of galaxy evolution. 

In § [2] we outline a simple semi-empirical way of estimat- 
ing the "intrinsic" BH mass/Eddington ratio distribution at a given 
AGN luminosity; we further demonstrate how dilution and radia- 
tive inefficiency can change this distribution. We also compare to 
observations in different wavebands and discuss how the observa- 
tional inferences relate to theoretical models (taking into account 
selection effects). In § [3] we examine how dilution and radiative 
inefficiency influence estimates of the "obscured" or Type 2 AGN 



fraction. In §|4]we investigate how the observed clustering of AGN 
as a function of wavelength, luminosity, and redshift can depend 
on the sample selection criteria, and in §[5] we examine how these 
effects can change the resulting host morphology distributions at 
a given AGN luminosity. Finally, in §[6] we summarize our results 
and outline observational tests that can resolve the degeneracies 
highlighted in this paper. 

For ease of comparison, we convert all observations to bolo- 
met ric luminosities given t he approp riate bolometric corr ections 
from lHopkins et al.l J2007bl) (see al s olRichards et alj|2006l) base d 
on the observations in lElvis et alj dl994l); iGeorge et al.l ( 1998); 



Vanden Berk et alj {200 1 ) ; Perol a et al.ld2002h ; lTelfer et alj (:2002); 
Ueda et alj d2003h : Ivignali et al] d2003l) . We ado pt a Q M = 0.3 , 
Qa =0.1, Ho =70kms"'Mpc"' cosmology and a lSalpeteildl955t) 
stellar initial mass function (IMF), and normalize all observations 
and models appropriately (note that this generally affects only the 
exact normalization of quantities here, not the qualitative conclu- 
sions). All magnitudes are in the Vega system. 



2 EDDINGTON RATIO DISTRIBUTIONS AND AGN 
SELECTION 

In order to construct plausible estimates for the effects of dilution, 
we begin with a model for the AGN Eddington ratio d istribution. 
We follow the methodology of Hopki ns et alj (2006a): given the 
observed AGN luminosity function, a model for the average AGN 
lightcurve (or lifetime) uniquely translates to an Eddington ratio 
distribution as a function of luminosity and redshift. For example, 
if lightcurves were step-functions with some "on" and "off" m, the 
Eddington ratio distribution would simply be the corresponding 5- 
functions. Here, we adopt a simple parameterization of the life- 
time/light curve motivated by hydrodyn amic simulations with BH 
feedback (Hopk ins et alj|200 5a. 2006b). The lightcurves, Edding- 
ton ratio distributions, a nd detailed comparison of bo th with obser- 
vations are presented in Hopkins & Hernquist| d2009h ; a number of 
other details and observational tests of the model (comparison with 
e.g. host mass distributions, clusteri ng, and other properties) are 
presen ted in a series of papers (e.g. iHopkins et al]|2005bl . 2006a. 
2008b). We refer the interested reader to these works for more de- 
tails, but briefly summarize the methodology here. 

We begin with the observed (redshift-dependent) bolometric 
quasar l uminosity function, he re from th e compilation of ob serva- 
tions in dHopkins et al.|[2007bl) (see also IShankar et alj|2009l and 
references therein). We parameterize the AGN lifetime (the statis- 
tical representation of the lightcurve; namely, the differential time 
spent per unit luminosity, for a BH of a given mass) with a simple 
Schechter-function fit: df/dlogL ~ to (L/Lo)~ a exp (— L/Lq). The 
normalization fo ~ 10 8 yr is the lifetime of hi gh-Eddington r atio ac- 
tivity and is comparable to the Salpeter time dSalpeteJl964l) for the 
e-folding of a BH at the Eddington limit; this estimate i s also rea- 
sonably consistent with observations (e.g.. lMartinill2004l and refer- 
ences therein). The "cutoff" at Lo ~ 0.4LEdd is the expected cutoff 
near the Eddington rate for the given BH mass, and a(Mmi) ~ 0- 5 
is a weak mass-dependent scaling (see Hopkins & He rnquistl2009h . 

The formalism used here can represent a range of realis- 
tic AGN light curve shapes, including those with an exponen- 
tial rise/fall, those with a range of variability power spectra, 
and those with a late-time power-law like de cay as predicted in 
a variety of self-regulated growth scenarios dMenci et alj 120031: 
i Granato et alj |2004|; ISazonov et alj |2004 IHopkins & Hemquistl 
boodlCiotti & Ostriker!l2007h . Adjusting the model parameters to 
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Figure 1. Top: Distribution of Eddington ratios (normalized dW/dlog A, where A = Lbol/^Edd) at a given AGN bolometric luminosity (as labeled), at z ~ — 1, 
derived from the models described in Left: "True" distribution - i.e. that in a "complete" sample that can recover all AGN regardless of radiative efficiency 
and/or dilution (appropriate for deep X-ray and host-galaxy selected deep narrow-line surveys). Middle: Distribution if only "radiatively efficient" AGN are 
detected, where we assume that thin disks disappear below m = 0.01. This is appropriate for surveys relying on broad lines and/or optical/IR continuum 
features (however deep). Right: Distribution if only AGN with optical/IR luminosities greater than ~ 10% of the host galaxy in the same band are selected. 
This is required for any color/continuum selection in these wavelengths, but also applies to most searches for broad-line signatures in complete galaxy 
samples (especially if the width of those lines is desired). Note the increasingly narrow range in the observed Eddington ratio distribution - and decreasing 
dependence on L — as these selection effects are included. This is because low-A AGN are potentially radiatively inefficient, and relatively less luminous 
relative to their hosts. Bottom: The observed Eddington ratio distributions for observed samples of AGN at a gi ven AGN bolometric luminosity (based on 
bolometric correction from the X-ray or IR luminosities). Left: X-Ray selected AGN from the observations in Hickox et al ] 120091) , where the effects of 
radiative efficiency and dilution are weak. This agrees well with the prediction for a "complete" sample at each Lb i- Center: Mid-IR sele cted AGN, where 
radiatively inefficient AGN will be under-luminous (and dilution effects enter as well, albeit to a lesser extent), also from lHickox et alj 120091) . The distribution 
is less luminosity-dependent, in r easonable agreement with predictions for such a sample. Right: Optically selected broad-line AGN from Kollmeier et al. 
(2006, atL bol > 10 44 ergs-') and lGreene & Hol l2007t at Lboi < 10 44 ergs '); for these samples radiatively inefficient AGN will not appear and AGN with 
luminosities much less than their hosts will also be excluded. Here there is almost no dependence of Eddington ratio on luminosity, with all sources having 
high Lbol/^Edd ~ 0.1 — 1, as predicted for samples with similar observational selection criteria. 



cover this broad range of possibilities, we find that our qualita- 
tive conclusions are relatively robust. Given this, we use the values 
quoted above based on hydrodynamic simulations. Likewise, vary- 
ing the choice of the bolometric QLF adopted by using just a subset 
of the observations compiled in Hopk ins et alj (12007b) makes little 
difference. 

Together, the bolometric luminosity function and prescription 
for AGN lifetimes specify the Eddington ratio distribution: the 
number at a given luminosity is known, and since the probabil- 
ity of a given system being at a given luminosity is trivially re- 
lated to the shape of the lightcurve, one can de-convolve to deter- 
mine the underlying BH mass/Eddington ratio distribution. Once 
this is specified, the most important aspect of the model is spec- 
ified. Given the observed Mbh — Mb u i B e relation (with normaliza- 
tion and scatter taken from lHaring & Rixl 12004). we can convert 



this to a host bulge mass distribution. Likewise, given the ob- 
served correlation between galaxy mass and bulge-to-total mass 
ratio (B/T), and its scatter (here from iBalcells et all 120071) . we 
can convert this to a distribution in total host galaxy mass and/or 
host morphology. Given some template galaxy SED or, more sim- 
ply, an average stellar mass-to-light ratio as a function of galaxy 
mas s (and redshift), well-determined in a number of studies (see 
e.g. iBell et alj|2003l; IConrov & Wechslerll2009l: IPbert et al.|l2009l : 
Ivan Dokkum et alj |2004 iKassin et alj |2007| ; iMoster et alj 120091) . 
we trivially convert the host mass to a host luminosity distribu- 
tion. Note that directly using the ob served correlation between BH 
mass and total host luminosity (e.g.jKormendv & Richstondl 19951) 
- accounting for its larger scatter - yields a nearly identical result. 

In Hopkins & Hernquistl (2009), we compare the distributions 
of BH and host properties inferred via the above procedure with 
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the distribution of Eddington ratios and BH and host masses mea- 
sured in deep narrow line samples at z = and inferred at z = 1 
(from iHeckman et alJl2004l : iHopkins et al]|2006cl : lYu et al.ll20ol 
iMerloni & Heinzll2008h . The agreement is good over the observed 
range from Mbh ~ 10 6 — 10 9 Mq. Likewise, comparing with the 
host galaxy mass , luminosity, and color distribu tions in these low- 
redshift samples (Kauffmann & Heckman 2008) and smaller high- 
redshift samples (e.g. Peng et al. 2006) yields reasonable agree- 
ment. So these simple analytic models at least provide a good ap- 
proximation to the actual observed distributions in the local Uni- 
verse. 

Given the above models, we take into account the effects of di- 
lution and/or radiatively inefficient accretion by considering three 
simplified models for how the "true" distribution might be ob- 
served: 

• (1) Complete: All sources at the given bolometric luminos- 
ity are included, regardless of Eddington ratio or host properties. 
This might be representative of e.g. deep X-ray surveys (modulo 
Compton-thick sources) or narrow-line optical surveys. 

• (2) Brighter than Host/Dilution-Limited: Given the observed 
depen dence of galaxy m ass-to-light ratios on stellar mass in various 
bands iBell et alj|2003h . we convert our host bulge plus disk mass 
distribution to a host luminosity distribution as noted above, and 
compare it to the observed AGN luminosity in the same band (with 
bolometric corrections from lHopkins et al]2 007b). We include ob- 
jects only above some threshold in the observed band; specifically, 
those sources for which the 6-band luminosity of the AGN is larger 
than 10% of its host luminosity in the same band. This is appro- 
priate for AGN samples selected on the basis of their broad-band 
IR-UV colors. 

• (3) Radiatively Efficient only: We include sources with Ed- 
dington ratios above some critical threshold, m > m C rit, where we 
take ra cr i t = 0.01 (see §[T). If there is indeed a change in accretion 
physics around this accretion rate, then such a cutoff is appropriate 
for samples that select AGN based on broad-band luminosities or 
colors in the UV, optical, or near infrared, or surveys that depend on 
identifying Type 1 signatures such as the presence of broad emis- 
sion lines. 

Figure [T] shows the inferred distribution of Eddington ratios at 
various bolometric luminosities for the three different sample cuts 
above[3 At a given luminosity (unlike at a given BH mass), the Ed- 
dington ratio distribution is roughly a lognormal distribu tion, albeit 
with some non-negligible skewness (as inferred by e.g. iFine et all 
2008). This is because, at sufficiently low m, arbitrarily high-MsH 
BHs would be implied (at fixed L), but the number of such sys- 
tems is vanishingly small. Figure Q] shows that the Eddington ratio 
distribution can be dramatically affected by different observational 
limits, becoming much narrower and more weakly dependent on 
luminosity when luminosity/Eddington ratio cuts are applied. 

Figure [2] shows the median and la dispersions in Eddington 
ratio as a function of luminosity at z — 0.5, based on approximating 

3 In order to be consistent with typical observational procedures, the "bolo- 
metric luminosity" plotted here is technically the bolometric luminosity 
that would be inferred based on the actual luminos ity in some band (given 
the bolometric corrections from Hopkins et al. 2007b, used throughout). If 
there is no qualitative change in SED shape with m, this is the true bolo- 
metric luminosity, but in the "radiatively efficient" sample cut above, this 
represents instead the inferred luminosity assuming, as is typically done, a 
standard high m SED, even if the system is in fact radiatively inefficient and 
thus possesses an intrinsically different SED. 
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Figure 2. Median Eddington ratio A and 1 a dispersion in Eddington ratio 
at a given bolometric AGN luminosity (as in Figure [TJ. We compare the 
predictions for the various selection criteria in Figure [T] with the observed 
distributions inferred from the distribution of X-ray to host luminosities 
in opt ically obscured AGN jFioreetal.1200 3; Hasinger 2008i; lHickox et alj 
black diamonds) (translating host luminosity to average BH mass), 
and fr om the distrib ution fitted to broad-lin e optical sample s inlFine et al 
1 2008) from the 2dF. lKollmeier et al J {moth from AGES, and lGreene & Hc 
12007) using the optical virial (line-width) BH mass estimators. At high-L, 
all the selection criteria yield similar results in good agreement with the 
observations - almost all systems at these L are radiatively efficient and 
bright relative to their hosts. At low-L, the X-ray and broad-line results 
disagree - the low mean A and wide dispersion in A near L ~ 10 43 ' 5 — 
10 44 ergs~' (L2-l0keV ~ 0.3 — 1 X 10 43 ergs~') for the X-ray sample is 
consistent with the predictions for a complete sampl e, while the large A ~ 
0.1 and small dispersion in the broad-line sample of Greene & Ho] j2007ri 
agrees with the predictions for a restricted (dilution or efficiency-limited) 
sample. 

the predicted distributions as lognormal (we use the IPV width to 
prevent bias from outliers or skewness). We compare the predic- 
tions _wth_several_observational estimates. 

IFine et al] d2008l) consider the Type 1 quasar population near 
z ~ 1 in the 2dF survey, and estimate the distribution of BH masses 
in narrow bins of luminosity relying on the commonly adopted 
virial BH mass estimators (see e.g. IVestergaard & Petersonll2006l 
and references therein). This allows them to consider the distri- 
bution of BH masses via this proxy as a function of luminosity 
down to near Seyfert luminosities. Because the width of the dis- 
tribution can be determined without relying on the (systematically 
still uncertain) absolute normalizations of these calibrators, the au- 
thors decline to estimate absolute Eddington ratios (although a 
rough estimate sugges t s the y lie between ~ 0.1 — 1, as predicted 
here). iKollmeier etafl ( 2006) use the same technique over a nar- 
row luminosity range of Type 1 AGN in the AGES survey. At yet 
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lower-luminosities, I Greene & Hoi d2007l) study a sample of local 
broad-line AGN from the SDSS down to extremely low luminosi- 
ties Lboi ~ 10 43 ergs' 1 (Lhc ~ 10 6 Lq). 

X-ray observations can study Type 2 objects at compara- 
bly low luminosities, where the optical light is dominated by the 
host galaxy and therefore a host galaxy stellar mass (and corre- 
sponding BH mass, ad opting the observed Mbh — relation from 
M arconi & H unt 2003) can be estimated. It is well-established that 
in this regime the optical luminosity of the galaxy is approxi- 
mately constant while the X-ray luminosity changes, implying that 
at lower-L the X-ray lumin osity function be c omes largely a se- 
quence in Eddington ratio (Fiore et al. 2003; Barger et al. 1 20051; 



La Franca et al.l 120051; ISilverman et al. 2005b a, 2008b]; lHasingerl 
2008tlHickox et alj|2009ri . " 



Figure [2] shows the results of the X-ray sample of iFiore et al.l 

d2003h . where the optical R-bsmd luminosity is converted to a stellar 
mass ba sed on the age and m ass-d ependent observe d mean M/L 
ratios in lBell & de Jond d2000h and lBell et all d2003l). We have re - 
calculated these comparisons using the samples of Hasinger (2008) 
and lHickox e"t~aH (2009) and obtain the same result; changing the 
assumed host M /L within uncertainties makes little difference. 

Figure [2] shows that the faint X-ray results agree with opti- 
cal results at high luminosities (as do the models). At low lumi- 
nosities, however, the X-ray results find that the median accre- 
tion rate decreases and the scatter increases, in agreement with our 
predictions for a "complete" sample^ Narrow-line AGN samples 
also find a similar tren d to that predicted for the "complete" sam- 
ple. lHao et all d2005al lbh find that host galaxy luminosity remains 
^constant as AGN lum inosity declines, down to very low luminosi- 
ties Lboi ~ 10 42 erg s~ 1 . Kauffm ann et alj d2003l) make this explicit, 
with Eddington ratios distributions moving to lower m < 0.01 and 
spanning an ~ 1 d ex range at the lowes t luminosities. 

However, the lGreene & Hol d2007l) optical results are quite dif- 
ferent: they find a median Eddington ratio m ~ 0. 1 with a relatively 
small (factor ~ 2) scatter, similar to what is found in the optical 
broad-line samples at three orders of magnitude higher luminosity. 
At the same bolometric luminosity, the median in is much higher 
than that in X-ray selected samples and the spread is smaller. This 
result agrees well with the expectation from our selection effect- 
limited models: because only sources with clear optical broad-line 
signatures of a "quasar" are selected, low-Eddington ratio systems 
which are simply underluminous relative to their hosts (and may 
also be radiatively inefficient) are excluded. 

Figure [TJ directly compares the Eddington ratio distributions 
observed in different wavebands and the model predictions, for 
AGN spanning three o rders of magni t ude in bolometric luminosity. 
In the X-ray sample of lHickox et alj d2009l) . the observed distribu- 
tions agree well with the expectation from the "complete" model 
- in particular, the median Eddington ratio clearly shifts to lower 
values with decreasing luminosity in the manner predicted. On the 
other hand, in infrared or optically-selected samples, the results are 
similar to that seen in Figure [2] observed Eddington ratio distri- 
butions are much narrower and less luminosity-dependent[f] This 



4 We expect this to be independent of whether or not Compton-thick 
sources are missed in hard X-ray samples, since we are here assuming that 
whatever "true" obscuration is present is not host-property dependent. This 
may be uncertain, but if it were not so, it would only strengthen our conclu- 
sions. 

5 The observed infrared distributions are somewhat broader than the pre- 
dictions; most of this can be attributed to scatter in the bolometric correc- 
tions and BH mass estimators. The predicted distributions would also be 
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Figure 3. AGN SEDs: an unobscured, Type 1 QSO (from Hopkins et al. 
2007b), and that spectrum obscured by a column with Nu = 10 22 cm~ 2 
(with M ilky Way gas-to-dust ratios and an SMC-like reddening curve from 
|Pelll992l and a small ~ 5% optical/UV scattering c omponent); we also 
show the predicted spectrum of an optically thin RIAF I Ouat aert & Naravai] 
1999) and a model with an inner RIAF and outer (truncated) thin «-disk. 
Without detailed SEDs, the RIAF models are indistinguishable from obscu- 
ration in the optical, and yield power-law like X-ray spectra with hardnesses 
similar to obscured Type 1 objects. The luminosities are arbitrarily scaled 
to be comparable to the systems of interest in X-rays. 



is more pronounced for the optical sample than the infrared one, 
but this too could be anticipated. In the optical, dilution is more se- 
vere than in the infrared; however, radiative efficiency, if important, 
should influence the results in both these bands similarly. 



IMPLICATIONS FOR 
POPULATIONS 



OBSCURED" AND TYPE 2 



Given the significant differences in the Eddington ratios selected 
in different samples, we now study whether this will affect other 
inferences drawn from AGN samples. We first focus on the im- 
plications for "obscured" and "unobscured" populations. Ideally, 
AGN would be observationally classified as obscured on the basis 
of detailed multiwavelength SED estimates - in particular X-ray 
spectroscopy, which can definitively show an absorption cutoff (al- 
though even in this case, Compton-thick sources will not be iden- 
tified), and corresponding absorption spectra in the infrared. Al- 
though there are well-observed systems where the presence or ab- 
sence of a ctual dust/gas obscuration can be unquestionably estab- 
lished (e.g. lWevmann et al.ll98ll : lLaor et aljl997l : lzakamska et al.l 
2004), these observations are prohibitively time-consuming for the 
high redshift surveys that are used to try to statistically quantify ob- 
scured fractions as a function of e.g. luminosity and redshift. As a 
consequence, AGN in these samples are typically classified as "ob- 
scured" on the basis of simplified proxies: the absence of dominant 
AGN luminosity/colors (or a broad-line SED) in the optical, or X- 
ray hardness. If however, the AGN emission is not dominant in the 
optical, these might masquerade as "obscured" sources in observed 
samples. 

Figure [3] illustrates these issues. We first compare template 



broader if the transition to radiative inefficiency were more gradual than we 
model. 
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SEDs of a Type 1 quasar and a "true" obscured quasar. The ob- 
scured system is the product of the Type 1 template attenuated 
by a column of Nh = 10 22 cm -2 , with a Milky Way-like gas-to- 
du st ra t io, so lar metallicity, SMC-like dust reddening curves fitted 
in |Pe]fl992T). and photoelectric ab sorption in the X-ray follow- 



ing 



Morrison & McCammorJ dl983h : for details see lHopkins et al.l 



(2006a[). Differences in e.g. the gas-to-dust ratios and reddening 
curves of AGN will somewhat alter the quantitative details here, but 
our qualitative conclusions below are not sensitive to these choices. 

Dilution can easily hide an AGN in the optical or infrared; re- 
call (eq.O that at an Eddington ratio m ~ 0.01, even if the system 
is still radiatively efficient, the S-band luminosity of the AGN is 

< 10% of the host galaxy luminosity in a relatively old elliptical 
host galaxy (with high M */L and B/T)\ 'm a y oung spiral galaxy 
(with typical B/T < 0.2; IWeinzirl et all 120091) that might host a 
~ 10 6 - 10 7 MqBH, the fraction of the optical light drops to < 1%. 
At these ratios, the galaxy continuum dominates, and the system 
would be classified as Type 2 or obscured according to any color, 
continuum, or broad-band based photometric criteria. The "true" 
AGN SED is unobscured, but the optical/IR light is predominantly 
galaxy light (by definition, in this regime). Broad lines would in 
principle be detectable, but only with very deep observations - 
the peak broad line h eight in the optical SED , given the compos- 
ite quasar SED (see IVanden Berk et alj|200~ll) . would appear as a 

< 2% correction to the continuu m for the star-forming ca se above 
(comparing to model SEDs from lBruzual & Chariot 2003). This is 
prohibitively time-consuming for all but a few local, well-studied 
systems (even in the early type galaxy, reliable identification of a 
single broad line would require S/N ~ 30 spectroscopy; typical 
SDSS spectra, for comparison, are S/N ~ 6, leading to difficulty 
making broad-lin e identifications when the AGN fraction falls be- 
low ~ 30 - 50%; IVanden Berk et alj f2006). In practice, most ob- 
jects in the large samples used to classify systems as obscured or 
unobscured require that the AGN continuum dominate in the opti- 
cal in order to define a system as Type 1 (in detail, systems are often 
classified as Type 1 or 2 on the basis of whether a single galaxy or 
AGN template provides a better fit to the optical; in such a case, the 
requirement that the AGN dominate the light is explicit). 

Figure[3]further compares the Type 1 and true obscured SEDs 
with two template spectra for radiatively inefficient accretion mod- 
els. First, we consider a RIAF around a 5 x 10 8 Mq BH with a 
dimensionless accretion rate of m = 0.01. The dimensionless spec- 
trum vL v /Lem is a weak function of black hole mass at fixed m, 
and so can be trivially scaled to different luminosities (in Figure [3] 
we arbitrarily scale the SEDs to be comparable in X-ray luminos- 
ity, since we will be considering samples observed in fixed lumi- 
nosity bins). The spectral c alculations here are based on the mod- 
els of lOuataert & Naravarj a 1999T ). to which we refer the reader 
fo r details. The one im portant difference is that we use the results 
of ISharma etal] J2007I) for electron heating in RIAFs. One uncer- 
tainty in the modeling is that RIAFs have strong outflows that may 
carry away a significant fractio n of the mass supplied at large radii 
dBlandford & Begelman|[l999l) ; for the models of interest here, the 
emission is largely produced at small radii close to the black hole 
horizon, and so the spectrum is not strongly sensitive to the pres- 
ence of outflows so long as the accretion rates quoted here are inter- 
preted as m at small radii close to the black hole's horizon; equiv- 
alently, we can consider systems at fixed luminosity regardless of 
the large-scale accretion rates producing that luminosity. In short, 
this affects only the ratio of the mass accretion rate to luminosity; 
since we will consider objects at fixed luminosity, this makes no 
difference. 



At the accretion rates relevant to this paper, theoretical mod- 
els and observati ons of both X-ray b inaries (e.g. lEsin et aill 19971) 
and AGN (e.g. lOuataert et al.ll 19991 ) suggest that the inner RIAF 
is likely surrounded by an outer thin disk. The latter produces 
a non-negligible optical-IR flux, although much less than a stan- 
dard radiatively efficient disk would produce if it extended close 
to the black hole. In addition, the disk photons inverse Comp- 
ton cool the inner RIAF, softening the X-ray spectrum. This is 
shown explicitly in Figure [3] where we also plot the spectrum 
for a m odel in which the transition from thin disk (modeled as a 
standard Sh akura & Sunvae4l973l a-disk) to RIAF occurs at 100 
Schwarzschild radii. 

In either case, it is straightforward to see how the RIAF spec- 
tra would be observationally equivalent to an obscured AGN in the 
optical. Interestingly (as we discuss further below), the RIAF spec- 
tra are also hard power-laws in the X-rays: their hardness can be 
similar to a gas-obscured intrinsic Type 1 AGN spectrum. 

Given our models for the Eddington ratio distributions of 
AGN, we can estimate the fraction of systems that would erro- 
neously be classified as optically "obscured" because of either di- 
lution or, more speculatively, a transition to radiatively inefficient 
accretion. This is simply the fraction, relative to a "complete" sam- 
ple, that meet our previous criteria (2) or (3), respectively. Figure|4] 
shows the resulting "obscured" fraction as a function of luminosity 
and Figure [5] shows the "obscured" fraction as a function of red- 
shift. Because at lower luminosities the fraction of low-Eddington 
ratio systems is larger (see § [TJ, this generically leads to a de- 
creasing "obscured" fraction with increasing luminosity. Moreover, 
because the shape and break location of the luminosity function 
evolve with redshift (implying that the distribution of Eddington ra- 
tios shifts accordingly in high-mass BH s; recall here we ado pt the 
QLF evolution from the compilation in iHopkins et afll2007bl but 
other observations yield similar results), at fixed luminosity there 
can be a weaker but still non-zero dependence of the "obscured" 
fraction on redshift (Fig. [5}. Specifically, it is well-established that 
more massive BHs are preferentially active at higher redshifts; as 
such, it is possible that a population at fixed (sub-L*) luminosity 
includes a larger fraction of high-mass BHs at low Eddington ratio, 
leading to larger (fractional) dilution. In addition, for the same host 
galaxy mass, star formation rates (and therefore stellar population 
luminosities) increase with redshift. 

The trends we find are similar to the in creasingly well 
established ob s erved t rends in vari o us sam ples ( Hill et al .lll99i 
I Simpson etal] 1 19991; Iwillott et al] l200Ct Simpson & RawlingL 
200d; ISteffen et alj 120031; lUeda et alj 120031; iGrimes et alj |2004 
i Hasinged 120041: ISazonov & Revnivtsevl |2004 iBarger & Cowiel 
20051; ISimpsoiJ |2005t lHao et alj l2005bl ; IgUU et alj l2007t 
iHickox et alj 120071; lHasingerl |2008|) . We compare the obser- 
vations from iHasingen J2008T) . a hard X-ray selected sample, to 
our simple estimates in Figures [4] & [5] Although some of the 
observationally identified "obscured" systems are almost certainly 
obscured in the real, traditional sense (especially at the highest 
luminosities), our estimates suggest that a reasonable fraction 
of the very high inferred obscured fractions (approaching unity) 
may owe to these selection effects, rather than to actual dust/gas 
obscuration. 

In detail, even with the strong assumptions in Figures [4] & 
[5] that all objects below Lb, agn/ Lb. host < 0.1 (dilution-limited) 
or below m = 0.01 (radiative efficiency-limited) are classified as 
"obscured," we do not reproduce the full observationally inferred 
obscured fraction. Including a constant (luminosity and redshift- 
independent) "true" obscured fraction of 0.2, our inferred "ob- 
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log(Lboi) [ergs 1 ] log(L bol ) [ergs 1 ] 

Figure 4. Predicted fraction of AGN at a given bolometric luminosity that would be classified as "obscured" in optical samples. Top Left: Fraction of objects 
at a given bolometric luminosity that are radiatively inefficient and thus would be absent in opti cal sample s. Solid and dashed lines correspond to the allowed 
range of Eddington ratio distributio ns der i ved fr om models and observations (see Hopkins & Hernquist 2009). We compare with the observed "obscured" 
fraction from the observations in Hasinger 1 2008), where (X-ray selected) objects are labeled "obscured" based on the lack of dominant broad-line AGN in the 
optical. Top Right: Same prediction, but with a constant (luminosity-independent) fraction = 0.2 of genuinely obscured sources added at all L. The combination 
of a true 20 % obscured fraction and radiatively inefficient accretion at low L can approximately reproduce the observationally inferred "obscured" sources. 
Bottom Left: As top left, but showing the predicted fraction of systems which will drop out in the optical from a dilution-limited sample (or samples where 
Type 1/2 classification is based on optical/UV/IR colors, continuum, presence of broad-lines, or goodness-of-fit of template AGN versus galaxy SEDs in 
the optical/IR). Bottom Right: As bottom left, with a constant fraction = 0.2 of genuine obscured sources. The observed trend of "obscured" fractions with 
luminosity can be explained primarily by these Eddington ratio effects, with a true obscured fraction ~ 0.2 is that is much less dependent on luminosity than 
a direct interpretation of the data would imply. 

soured" population is similar to observational estimates as a func- 
tion of both luminosity and redshift (right-hand panels in Figs.|4]& 
PS}- Thus, although some genuine obscuration is required, it might 
be as low as ~ 20% of the population, significantly less than what 
is typically inferred in low-luminosity AGN samples; more impor- 
tant, the majority of the inferred luminosity and redshift depen- 
dence in observed samples can be accounted for by the effects of a 
non-trivial Eddington ratio distribution and the inevitable selection 
effects in optical continuum or broad-line samples. 

Note that we are not arguing that the obscured fraction has to 
be as low as ~ 20%. Given the uncertainties (error bars in the obser- 
vations, range of allowed Eddington ratio distributions and system- 
atic uncertainty in our simple model of host properties), it could be 
somewhat (factor ~ 2) higher (as suggested by other independent 
constrai nts from e.g. infrare d emission or X-ray background fitting; 
see e.g. iTreister et alj|2008h . What we show, however, is that sim- 
ple effects could, within reasonable parameter space, account for 
up to ~ 80% of the claimed effects, in particular, the strong de- 
pendence on luminosity and redshift. In fact, direct observations 



that do attempt to more rigorously distinguish obscured and un- 
obscured sources in "complete" samples (e.g. narrow-line studies) 
favor a somewhat higher "true" obscured fraction ~ 30 — 40%, but 
wit h a similar lack o f luminosity and/or redshift dependence (see 
e.g. lHao et all2005j|bl and references therein). 

Qualitatively similar observational results to those shown in 
Figure [4] on the fraction of obscured systems have been obtained 
by using X-ray hardness (rather than optical SEDs) as a proxy for 
obscuration; i.e, systems with sufficiently hard X-ray spectra are 
classified as obscured (although it should be noted, the statistics 
in these studies are much more limited). If real, this is difficult 
to explain with simple host galaxy dilution. However, in Figure fj] 
we noted that the expected RIAF spectra are similar in their X-ray 
hardness to a genuine obscured (intrinsically Type 1) SED. We con- 
sider this in more detail in the top panel of Figure [6] which shows 
the X-ray SEDs for Type 1, genuinely Type 2 dust/gas-obscured, 
pure RIAF, and RIAF+outer thin disk models.The RIAF SEDs are 
hard power laws; they do not have an absorption cutoff. But, if we 
ignore that cutoff, only seen at < 1 keV rest-frame, the RIAF SEDs 
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Figure 5. As in Figure [4] but plotting the predicted fraction that would 
appear "obscured" in optical samples at a given L (as labeled) as a function 
of redshift. Also shown are the corresponding observational estimates from 



Hasinser 1 2008) for luminosities in the range L ~ 10 — 10 erg s 



are sufficiently hard that they are indistinguishable (with typical 
signal-to-noise) from the SED of the true obscured AGN; both are 
significantly harder than the SED of an unobscured AGN out to 
> lOkeV. 

To quantify the extent to which our theoretically calculated 
RIAF spectra could mimic obscuration, Figure [6] also shows the 
hardness ratio (HR) as a function of galaxy redshift for these tem- 
plate SEDs, with hardness ratio defined in terms of the observed- 
frame hard (2-10 keV) and soft (0.5-2 keV) band counts, HR= 
(H-S)/(H + S)E 

Larger absolute values of HR correspond to 
"harder" SEDs and, it is usually believed, a more obscured system. 

The RIAF model, being a near power-law, predicts a con- 
stant, relatively hard HR. We compare this to the HR derived from 
the Type 1 template shown, obscured by different columns Nh = 
10 22 , 10 23 , 10 24 ergs~'. The templ ate shown corresponds to an in- 
trinsi c photon index T = 1.9 (see [Elvis et al J 1 1 9941: iGeorge et al.l 
1 19981 : IPerola et alj 120021 : lUeda et alj 120031 : Ivignali et alj l2003h . 
with a reflection component generated b y the PEXRAV m odel 
dMagdziarz &Zdziarskilll995h following lUeda et all J2003h . We 
also show the results for varying T, within the range suggested by 
observations, V — 1.8 — 2.0 (dotted lines in the middle panel of Fig. 
1. 

At z > 1, the strong absorption cutoff in true obscured sys- 
tems is redshifted out of the soft band, and the hardness ratios from 
obscured systems are similar to the power-law spectra expected in 
RIAF models (middle panel; Fig. [6}. We therefore consider what 
column densities would be inferred, if the true spectrum were the 
RIAF model observed at some redshift, but the observed hardness 



6 Note that the observed hardness ratio, as defined in terms of counts, de- 
pends on the exact waveband and instrument response function. For the sake 
of generality, what we quote here can be thought of as an "intrinsic" HR for 
an instrument with perfect response and sensitivity, with fixed bands defined 
from 0.5 — 2.0 keV and 2.0 — 10.0 keV; as such, it is not directly compara- 
ble to hardness ratios in the literature for real instruments and wavebands. 
Because we are interested in comparing two input SEDs (the obscured AGN 
and RIAF), however, so long as we treat them identically it makes little dif- 
ference to our conclusions whether we adopt this definition or model a more 
complex instrumental response function. There is, however, some weak sen- 
sitivity to the wavebands in which hardness ratios are defined, since these 
will sample various portions of the absorption cutoff at different redshifts. 
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Figure 6. Top: Comparison of template AGN SEDs in the X-ray; we con- 
side Type 1, obscured, pure RIAF, and RIAF+outer disk models. The RIAFs 
predict hard X-ray SEDs similar to an obscured system, but without a true 
absorption cutoff. Middle: Observed hardness ratio (HR) from the templates 
in the top panel as a function of redshift. The Type 1 template is obscured 
by different columns Nh = 10 22 , 10 23 , 10 24 ergs _1 for input photon index 
T = 1.9 (solid lines) and T varying from 1.8 - 2.0 (dotted lines). At z > 1, 
the strong absorption cutoff in true obscured systems is redshifted out of 
the soft band, and the hardness ratios from obscured systems are similar 
to the power-law spectra expected in RIAF models. Bottom: Inverting the 
middle plot: Given the template RIAF hardness, we plot the column density 
Nh that would be inferred by matching this hardness to an obscured Type- 1 
SED. We assume observatio nal errors of ~ .1 in HR, comparable to deep 
observed X-ray samples (see Hasinger 2008). 
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Figure 7. Model RIAF spectra at Eddington ratios m = 
0.025, 0.01, 0.005, 0.0025, 0.001 (solid; black, red, orange, green, 
blue, respectively, from top to bottom in the X-rays). Also shown are 
models with m = 0.01 in which there is an outer thin disk that transitions to 
the inner RIAF at a radius of 10 (dotted) and 100 (dashed) Schwarzschild 
radii. For the simple comparisons made in this paper ("obscured" or 
"unobscured"), the differences between the different m have little impact; 
the optical SEDs are always strongly suppressed, and the X-ray SEDs 
are power-law like (over ~ 0.5 — lOOkeV) with similar hardness for the 
accretion rates of interest, down to very low-L (L DO i <C 10 42 erg s~ 1 ) where 
the spectra become softer. 



ratio were modeled as the template Type 1 SED plus obscuration 
(allowing for typical errors in hardness ratio AHR ~ 0.1, typical 
of the quality in deep X-ray samples). This is th e common method 
for inferring obscuration in X-ray samples (e.g. lUeda et al.1120031 ; 
lLa Franca et al.l2005l : lBarger & Cowiel2005l : lGilli et al.l2007l) ; ap- 
plied to a hard power-law like spectrum expected in radiatively in- 
efficient states, it can lead to a significant false obscuration sig- 
nal, with inferred column densities Nh ~ 10 22 — 10 23 cm" 2 (bottom 
panel of Fig. [6j. At these redshifts, observations are sampling the 
regime where the RIAF SEDs are significantly harder than that of 
an unobscured AGN, but missing the regime where the strong ab- 
sorption cutoff distinguishing true absorption can be sampled. 

Figure|7]shows how the RIAF spectrum varies as a function of 
accretion rate. As rh decreases, the predicted X-ray spectra become 
softer (a transi tion that has also been suggested from observed 
AGN ; see e.g. iDai et aljpiool ; ISaez et aljHobljl : IConstantin et~aH 
2009); however, they are still harder than an unobscured AGN 
for all m > 0.001 (below which the luminosities are sufficiently 
low that they contribute negligibly to the different fixed-luminosity 
samples considered here). It therefore makes little difference 
whether we simply pick one of these Eddington ratios as "represen- 
tative" of the RIAF SED, or include the full m dependence in the 
predictions here. Changing other parameters in the model shown 
in Figure|7] such as the disk viscosity, the magnetic energy density 
(relevant for the predicted synchrotron emission), and mass loading 
of outflows does not significantly change the spectral shape at fixed 
rix. The strength of outflows is important for the normalization of 
the SED, i.e., for the ratio of the luminosity to the (large-scale) ac- 
cretion rate, but at fixed luminosity makes little difference. In other 
words, although some of the details of RIAF accretion flows remain 



uncertain, the vast majority of these issues do not affect our quali- 
tative conclusions. RIAF SEDs are generally X-ray hard; this must 
be true, since they are designed to describe the observed "low/hard" 
(low-luminosity, hard X-ray spectrum) state in X-ray binaries. The 
dominant uncertainty for our purposes lies in the nature of the tran- 
sition to the RIAF state. If this is more gradual than what we have 
assumed here - i.e. if a significant thin disk survives all the way 
down to small radii for Eddington ratios m < 0.01, then the X-ray 
SED could be significantly softer (owing to inverse Compton cool- 
ing of the inner thick disk photons), although it must survive to very 
small radii indeed (only being an important effect if the transition 
radius is < 10 Schwarzschild radii). Of course, this is essentially 
equivalent to the uncertainty re garding wh ether such a transition 
occurs at all (see, for example. lMaodl2007r) . 

It is also important to note, as is clear in Figure [6] that the 
RIAF SEDs do not become extremely hard. It is difficult, with 
such models and no true obscuration, to mimic e.g. Compton thick 
columns with Nh > 10 23 cm" 2 . Ultimately, more detailed observa- 
tional studies will allow direct comparison of the hardness ratio 
distribution at each redshift, which could include components that 
cannot be mimicked by RIAF models. 

We now estimate the fraction of AGN that might be classified 
as "obscured" on the basis of observed X-ray hardness, as a con- 
sequence of having an intrinsic RIAF-type SED. This is shown in 
Figure[8]as a function of luminosity and redshift. The model is the 
same as in our predictions for optical samples, but we assume that 
all sources with m < 0.01 have a thick disk spectrum as in Figure[3] 
and calculate the hardness ratio distribution and implied column 
density di stribution as in F igure [6] As in the observational litera- 
ture (e.g., iGilli et alj|2007l and references therein), we classify as 
"obscured" all sources with an inferred column Nh > 10 22 ergs _1 . 
The radiatively efficient {m > 0.01) sources are assumed to be ei- 
ther uniformly unobscured ("true /obscured = 0" case; left panel) or 
to have a constant fraction /obscured = 0.2 of sources genuinely ob- 
scured (right panel). The resulting trends in obscured fraction as 
a function of luminosity and redshift are similar to those in Fig- 
ures |4]& [5] (the optical samples). By taking into account the pre- 
dicted change in AGN SEDs below m ~ 0.01, we thus conclude 
that a significant fraction of the inferred dependence of obscuration 
on luminosity/redshift, both determined from optical observations 
and X-ray hardness cuts, could owe to the fraction of systems in a 
radiatively inefficient state, rather than from true obscuration. 



4 AGN CLUSTERING VERSUS SELECTION AND 
FUELING 

By affecting the distribution of Eddington ratios detected in a given 
sample, dilution and radiative efficiency selection effects can also 
affect the BH masses to which the sample is sensitive, and thus 
stellar/host masses. Given the strong dependence of clustering on 
stellar/halo mass, the overall clustering amplitude measured in such 
a sample can be different than that inferred from a sample including 
all active systems. 

Another important effect for observed clustering is related to 
the significant time systems can remain active after their initial 
peak of activity; if, after a near-Eddington peak in activity, the 
AGN accretion rate declines in an approximate power-law fashion 
tx (t/to)^ 13 , then for a system at much lower Eddington ratio m, 
it has been a time ~ tom~ [ ^ since the peak of activity. For char- 
acteristic to ~ 10 8 yr suggested by observational estimates of the 
AGN lifetime and f3 ~ 1 — 3 from a variety of different models (see 
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Figure 8. The fraction of systems that would be classified as "obscured" based on their observed X-ray hardness ratios (as shown in Figure [6] as a function 
of luminosity and redshift); we assume tha t for m < 0.01, all accretion is via a RIAF. The predictions are compared with the observed fraction of hard X-ray 
sources, typically classified as "obscured" lHasinger 2008). Because RIAFs have intrinsically hard spectra, X-ray hardness alone does not necessarily break 
the degeneracies between selection effects and true obscuration demonstrated in Figures I4l5l X-ray spectra with the sensitivity to measure the presence or 
absence of a true absorption cutoff are needed to definitively show obscuration. 



Hopki ns & Hemquistl200f3.l2009L and references therein), this time 
is still small compared to the Hubble time for high Eddington ratios 
m > 0.1. The system is therefore viewed with no significant evo- 
lution in its properties between the peak of activity and the time at 
which it is observed. However, at sufficiently low Eddington ratios 
m <C 0.01, the initial trigger may have occurred at a much earlier 
time > Gyr (comparable to a Hubble time), implying that systems 
observed at e.g. z = - having decayed for Gyr down to these 
Eddington ratios - were active/forming at much earlier times, and 
therefore are more biased than a system of similar mass first be- 
coming active at z = 0. Modeling this in detail requires a detailed 
and specific mo del for AGN lightcurv es an d triggering; we adopt 
the results fr om lMcOuinn et alj J200gh fsee lHopkins et al.ll200'M 
2008a,. l2009l for more details) to demonstrate the effects here; how- 
ever, many of the qualitative results are general. That model uses 
identical assumptions to the model described here in §[2] (the same 
AGN luminosity function, light curves, Eddington ratio and host 
property distributions), but places the mock AGN in halos in a cos- 
mological simulation according to the obs ervationally constra ined 
galaxy mass-halo mass relations (see e.g. iMoster et alj|2009l . and 



references therein), and so can predict clustering and other proper- 
ties. 

Figure [9] shows the bias as a function of luminosity taking 
into account the above effects and for our three different "obser- 
vational" samples. At the bright end, there is no difference between 
models and a reasonably steep dependence of clustering on lumi- 
nosity, because essentially all models agree that AGN must be near- 
Eddington at these luminosities, implying that the luminosity func- 
tion is a sequence in BH/host mass. At luminosities < L„, however, 
the fact that non-trivial lightcurves allow the population to be a mix 
of both high and low Eddington ratio systems rather than a strict 
sequence in mass makes the dependence of clustering amplitude 
on luminosity relatively weak, as is observed JPorciani et al. 1 20041; 
l Adelberger & Steidell Eool; IPorciani & Norberd 120061 ; ICoil etail 
120071 ; Ida Angela et alj|2008h (for more details regarding the mod- 
eling of clustering, see iLidz et al.| [2006). For "complete" samples, 
the clustering amplitude is predicted to increase at lower luminosi- 
ties, because of an increasing contribution of very massive systems 
that have decayed to very faint luminosities (low Eddington ra- 
tios). Indeed, observations of faint X-ray sources appear to support 
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Figure 9. Dependence of clustering on luminosity at z ~ for the different 
samples described in Figure[T] At L > L„, the samples are all similar, with 
mostly near-Eddington systems; thus L scales roughly with halo mass and 
there is a steep b(L). At lower L, the broad range of Eddington ratios leads 
to a weak dependence of b(L) and can cause differences in different samples 
(see lLidz et al.ll20061) . At very low L, the clustering in a complete sample 
can increase again (depending on the exact hghtcurve/lifetime model), as 
the luminosity function can be dominated by early-forming highly clustered 
BHs that have had time to decay to very low luminosities. These low m sys- 
tems would not appear in radiatively efficient or dilution-limited samples. 



this, wit h a surprisingly high cluster i ng amplitude in several cases 
(see e.g.lConstantin & VogelevfeOOoilPlionis et alj200&lciiiii et all 



l2009l : lHickox et al.l2009h - higher than observed in optical samples 
at the same redshift and similar or higher lumin osities dCroom et al.l 
I2005I: iMvers et al]2007l ; Ida Angela etai]|2008h . 

If, however, these sources are selected against in optical/IR 
samples, then the clustering amplitude remains flat or decreases at 
low luminosities, reflecting the fact that the most biased sources 
(which are very massive, old galaxies that with AGN at very low 
Eddington ratios) are excluded, leaving a sample dominated by 
lower-mass systems at intermediate/high Eddington ratios. This 
may explain why the clustering amplitudes in X-ray and optical 
samples appear to disagree at low luminosities. 

Figure [Tol shows the clustering amplitude expected as a func- 
tion of redshift for our standard three samples of AGN at two dif- 
ferent intervals in bolometric luminosity. At a relatively high lu- 
minosity (top panel), which corresponds roughly to ~ L, in the 
AG N luminosity function over the redshift range z ~ 0.5 — 2 (see 
e.g. iHopkins et al]|2007bl) , the different selection methods yield 
nearly identical clustering amplitudes, in agreement with a vari- 
ety of observations. At low luminosities, however (Lboi ~ 10 43 — 
10 44 ergs~', corresponding to Z^-iokeV ~ 10 42 — 10 43 ergs _1 ), the 
bias versus redshift depends on the selection criteria (bottom panel 
in Fig. UOt . In particular, at these luminosities, AGN selected in 
dilution-limited samples are preferentially at high Eddington ratio 
and therefore, smaller BHs in low-mass, weakly-clustered hosts. 
The bias as a function of redshift is somewhat sensitive to the ex- 
act model we adopt, but the selection effects are in any case clearly 
not negligible - they can lead to a factor of ~ 2 difference in bias, 
which corresponds to up to an order-of-magnitude difference in the 
implied host halo mass at these redshifts and bias values. This is 
qualitatively similar to what is observed in radio versus X-ra y ver- 
sus infrared-selected AGN samples (e.g.. lHickox et al. 2009). 
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Figure 10. Dependence of AGN clustering on redshift in a intermediate- 
high luminosity sample (~ L* quasars at these redshifts; top) and a 
low-luminosity sample (corresponding to L2-l0keV ~ 10 42 — 10 43 ergs _l 
or Mb ~ —17; bottom), taking into account the different selection ef- 
fects described in Figure [T] For the ~ L, sample, where most obser- 
vations have been made, the different selection effects do not signifi- 
cantly change the predicted bias, which is also consistent with the ob- 
served clust ering of quas a rs of s imilar luminosit ies from optically s e lected 
samples in Croom et al. ( 200^ red squares), iPorciani &Norberg (2006, 
green diamonds ), Myers et al I M20071 2006, cyan and blue circles), and 
Ida Angela eta?] (2008, violet stars). At very low luminosities, the effects 
described in Figure|9]and §|5]lead to significant differences between sam- 
ples; the higher bias in X-ray samples has been observed in a number of 
studies (see text). 



5 DIFFERENCES IN HOST MASSES AND 
MORPHOLOGIES 

At a given luminosity, different Eddington ratios imply different 
BH masses and thus different host masses. Given the strong corre- 
lation between galaxy morphology and stellar mass, this naturally 
leads to variations in host morphologies found in different samples. 
Recall, our model includes both the host total mass and bulge/disk 
mass distributions. For convenience, we therefore adopt the bulge- 
to-total ratio B/T as our morphological proxy of interest. 

Some caution is warranted here - making such a conversion 
implicitly assumes that whatever fuels objects is statistically ran- 
dom at each bulge mass/luminosity, and does not select out a spe- 
cific morphological sub-population at a given stellar mass. This is 
not necessarily true if e.g. mergers or disk instabilities are the dom- 
inant means of AGN fueling. Nevertheless, some samples suggest 
that, at least at low luminosities where fueling an AGN requires 
only a small gas supply, the host distribution spans a wide range 
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Figure 11. Model distribution of host galaxy stellar masses and B/T (bulge- 
to-total stellar mass ratio, our proxy for morphology) as a function of bolo- 
metric luminosity for the different selection samples considered in Figurefl] 
In a complete sample (X-ray or narrow line), even low-L AGN populations 
are dominated by relatively massive (M ~ IP 10 ' 5 - 10" Afo) early-type 
(Sa/SO/E) hosts, as is observed (e.g. Kauffmann et al. 2003). In a radiatively 
efficient or dilution limited sample, however, the requirement of moderate 
Eddington ratios at low L leads to host masses and morphologies that are 
more luminosity-dependent: lower-L AGN in these samples will be hosted 
by lower-mass, late-type (B/T < 0.3) galaxies. This can explain the ten- 
sion between the "conventional wisdom" that Seyferts live in disks (from 
optical-IR samples) and the contradictory results from recent narrow-line 
or X-ray samples that moderate-luminosity AGN have primarily early-type 
hosts. 
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Figure 12. Top: Median bulge-to-total ratio as a function of luminosity, 
for the complete (left) and dilution-limited (right) samples in the model. 
Dotted lines show the ± 1 a range in the population. Bottom Left: Corre- 
sponding observational estimates from the narrow-line SDSS AGN sample 
(Kauffmann et al. 2003), converting Loin t o bolomet ric and concentrations 
to B/T given the mean observed relation in Graham (2001) (error bars are 
the dispersion in the population, not errors in the median). Bottom Right: 
Same, based on the optically-selected PG quasar sample from Dunlop et al. 
1 2003). We show the results both for their sample with "robust" morpholo- 
gies (reliable bulge/disk decomposition; black circles) and their full sample 
(magenta squares). 



and does not show a strong preference for specific morphologi- 
cal types feahcall et a l. 1997; Malkan et al. 1998; Heckma netal] 
2004 iFalomo et alj|2004 ISanchez et al]|2004l ; iRigbv et alJbood : 
Nandra et alj|2007l) . At higher luminosities this is probably not the 
case, but as before, it is the low-luminosity regime where the Ed- 
dington ratio distribution is broad and has a significant effect. 

The predicted host stellar mass and morphology distributions 
are shown in Figure[TT] For the reasons discussed in §f4] a sample 
that includes objects of all types (including those that are radia- 
tively inefficient or diluted in the optical/IR) is dominated by rela- 
tively massive (~ L*) early-type galaxies at all luminosities. This 
conclusion has been reac hed observationally in e .g . narrow line 
and d eep X-ray samples: |K auffman n et alj ( l2003h . iKewlev et al.l 
d2006h and lHao et aU [2005b) show that host luminosity and mass 
are nearly independent of AGN luminosity, with early-type mor- 
phologies predominant (similar to the X-ray conclusions discussed 
in Figure[2]l. However, in a sample with an implicit Eddington ratio 
selection, only low-mass BHs are present in the sample at low lumi- 
nosities, which leads to a host sample dominated by less massive, 
more disk-dominated and star-forming host galaxies. 

This result can potentially resolve long-standing discrepancies 
between different observations of AGN host galaxies, in samples 
that should select for similar AGN bolometric luminosities. We il- 
lustrate this in Figure [T2l where we plot the median B/T as a func- 
tion of AGN luminosity for two model samples: our complete sam- 
ple and the dilution limited/optical broad-line sample (we do not 
plot host masses/luminosities because this information is already 
effectively presented in Figure(2] given Mbh oc Mb u i gc ; moreover the 
host morphology distribution is the primary focus of many stud- 
ies attempting to constrain AGN fueling processes). We compare 
with roughly comparable observational samples: the deep narrow- 



line sample of [Kauffmann et al. (2003) (which includes faint ob- 
jects with luminosities well below those of their hosts) and the 
broad-line PG quasa r sa mple (all of which h ave Lagn > i-hosi) from 
iDunlop et alj J2003h and lFlovd et all d2004l) H 

The "conventional wisdom," derived from early faint AGN 
samples that were dominated by broad-band optical/UV color 
selection, has been that Seyfer ts reside predominantly in disk- 
dominated, low-ma s s galaxies 1 Adamsl 19771 : [Smith et al, 1 19861 : 
lGranatoetai1ll993l ; iMcLeod & Riekdll995h : this is still true in 
modern samples that adopt similar selection methods, i.e., that re- 
quire a moderate UV/optical /IR luminosity from the AGN relative 
to the host galaxy (see e. g . iMartini & Poggelll999l : IPeletier et al.1 
Il999l : IDunlop et all 120031 : IZakamska et al]|200d and references 
therein). These inferences are consistent with our models for a sam- 
ple with an effective minimum Eddington ratio. However, suffi- 
ciently deep narrow-line and X-ray surveys ("complete" surveys) 
find instead that the AGN host population is dominated by massive 
early-type galaxies, at all luminosities. At a given (similar) bolo- 
metric luminosity, these surveys suggest that the AGN host fraction 
rises steeply as a function of bulge mass, in contrast to the results 
from Type 1-sensitive surveys. As Figure Q/2] shows, both results 
follow naturally once one takes into account a broad intrinsic dis- 
tribution of Eddington ratios and the observational selection effects 
highlighted throughout this paper. 



7 Note that the fitting procedures for B/T in these samples differ from the 
calibration in Balcells et al. ( 2007) that we use for our predictions - there- 
fore some systematic offsets are expected. The qualitative trends with mass 
are of primary interest here. 
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6 DISCUSSION AND CONCLUSIONS 

We have demonstrated that, given a non-trivial Eddington ratio dis- 
tribution, two simple effects will give rise to potentially important 
differences in the AGN populations selected by different observa- 
tional criteria. Even if the physics of accretion were independent 
of m, an AGN will inevitably be much less luminous than its host 
galaxy in a given band at sufficiently low Eddington ratio. Dilution 
is, of course, inevitable, and in most surveys AGN will be classified 
as "normal galaxies" if the AGN luminosity is less than that of the 
host. In the optical/IR, this happens at surprisingly high Eddington 
ratios, m ~ 0.01 — 0. 1, for most low-luminosity populations, which 
is why the effects of dilution can have such a significant impact 
on AGN surveys. Moreover, at low accretion rates m < 0.01, AGN 
may transition to a radiatively inefficient state, in which the optical 
and infrared emission decline rapidly, and eventually the broad-line 
region may disappear entirely (see footnote 2 in ^T]for more discus- 
sion of how the broad line emission disappears). Regardless of the 
specific spectral changes that may occur, the phenomenon of low- 
efficiency accretion is increasingly well-established and at low ac- 
cretion rates populations of AGN with little to no broad-line emis- 
sion appear dHoll 999; Mora n et all200d:lTrart2003l;lNicastro et al.1 
l2003l ; lBrightman & Na ndra 2008; lBianchi et alj2008l) . 

At high luminosities, observed systems are uniformly at rela- 
tively high Eddington ratios, and so it is unlikely that either dilu- 
tion or radiative inefficiency is a significant issue. At low luminosi- 
ties, however, the "true" Eddington ratio distribution of all AGN, as 
constrained by sufficiently deep narrow-line, X-ray, and radio sur- 
veys, is broad and extends to lower m with decreasing bolometric 
luminosity (roughly as A = Lbai/jjMd cx Z ^~' at sufficiently low 
luminosities; see, e.g., iFiore et al. 2003; Si lverman et all 12005a; 
Hicko x et al . 2009); this is consistent with some theoretical models 
of AGN lightcurve evolution ( Hopkins & Hernquist 2009). 

Given the selection effects implicit in the construction of a 
broad-line or Type 1 optical/mid-IR AGN sample, our results imply 
that dilution and radiative efficiency lead to these samples prefer- 
entially sampling the high-m end of the population, with rh > 0. 1 
and a relatively narrow range in m even at low luminosities. This is 
indeed observed (Figures l"H2l l. The strong impact of selection crite- 
ria is illustrated by the distinct Eddington ratio distributions recov- 
ered in X-ray samples (where dilution is weak, and radiatively effi- 
ciency ^ecUnesJessrapidly) as compared to optical samples (com- 
pare iGreenei&H^llQO^) . 

For the same bolometric luminosity, different Eddington ra- 
tios imply different BH masses and, correspondingly, host galaxy 
masses. Low-luminosity AGN hosts in sufficiently deep narrow 
line/X-ray/"complete" samples should cluster around a narrow 
range in host mass, pr eferentially massive (~ L«), early-type gala xy 
hosts, as observed jKauffmann et al.1 120031 ; Iffito et : all 12005 dlbl: 
IKewlev et alj 120061 ; Isilverman et alj l2008al : Iffickox et al J 120091) . 
However, low-luminosity AGN in dilution (or radiative efficiency)- 
limited samples will often be found in less massive hosts, which 
are correspondingly less bulge-dominated (Figures [l HI 2t . This is 
also observed, and can reconcile conflicting claims in the litera- 
ture regarding the hosts of Seyfert galaxies: samples based on op- 
tical/UV excess selection have traditionall y found that these ob- 
jects reside in less mass i ve late-type disks {B e Robertis et al .1 1 9981 : 
Martini & Poggej Il999l ; iKnapen et all l200d: iDunlop et al] 120031 ; 
Rigbv et al.l l2006l) while more recent X-ray or narrow line sur- 



fects also lead to differences in the expected clustering (Figure [9}: 
faint X-ray/narrow-line/"complete" AGN samples should exhibit 
higher clustering amplitude than corresponding optical samples. 
This is consistently observed when comparing X-ray and opti- 
cal/IR AGN clustering p r operties at low luminosities (see e.g . 
IPorciani & Norberd [20061 : iMvers eTail 120071 ; Iplionis et al.1 120081 
and references therein). 

In most surveys, AGN are identified as "obscured" or 
"unobscured" in the optical based on a simple classification, 
namely whether the optical spectrum is dominated by a broad- 
line AGN SED or galaxy light; we have shown, however, that 
both dilution and radiative inefficiency will lead to the clas- 
sification of an AGN as obscured (e.g., Fig. [3}. The natural 
expectation is that the fraction of sources which are optically 
dim for these reasons should increase with decreasing luminos- 
ity. It is in fact traditionally inferred that the intrinsically "ob- 
scured" (non Type 1) AGN fraction behaves in just this manner 



scurea (non type 1 ) AUJN traction behaves in just this mannei 
dffillet all 19961; I Willott et al.l200d:ISteffen et al j2003l;IUeda et al 



2003; Grimes et al 



20041; iHasingei 



2004; Sazonov & Revnivtsev 



2004l;fBarger & Cowiell2005l ; lGilli et al.ll2007l ; Iffickox et alJ^Oa 
lHasingerll2008l) . We find, however, that the fraction of AGN that 
are diluted (and may also be radiatively inefficient) scales in a sim- 
ilar manner to the observed optically "obscured" fraction, with both 
luminosity and redshift. The effect is not large enough to account 
for all optically dim objects (in particular at high L), but the frac- 
tion of "true" obscured sources could, in principle, be as low as 
~ 20 — 40%, with a much weaker dependence on both luminosity 
and redshift than that traditionally inferred (see Figures l4l8t . 

Using X-ray hardness as a measure of obscuration represents 
an improvement over the "optically dim" criterion in many ways 
(e.g., it should not be affected by simple dilution). However, the 
expected X-ray SEDs of AGN in radiatively inefficient states are 
relatively hard power-laws for m ~ 10~ 2 (Figure |7J- They are 
sufficiently hard, in fact, that the systems would be classified as 
obscured/hard sources by the standard observational criterion - 
which uses a simple hardness ratio to probe the presence or ab- 
sence of an absorption cutoff, given the lack of full X-ray spec- 
tra at the signal to noise ratios available for most high redshift 
sources. Convolving with the observed Eddington ratio distribu- 
tion, we find that RIAFs could explain a large fraction of the X- 
ray "obscured" sources at low luminosities (Lboi < 10 45 ergs~', or 
i-2-iokcV J; 3 x 10 43 ergs _1 ). In particular, just as with the optical 
samples, we find that the "true" obscured fraction could be as low 
as ~ 20 — 40%, with a much weaker dependence on luminosity 
and/or redshift than has been suggested by observations based on 
hardness ratio distributions. It is important to stress that we are not 
questioning the fact that local samples of Seyferts indicate a large 
population of obscured sources. Rather we are pointing out that the 
hard X-ray spectra detected in higher redshift samples may in many 
cases be intrinsic (having to do with accretion physics), rather than 
indicative of absorption. 

Consistent with the conclusions drawn here, detailed ob- 
servations have revealed that dilution accounts for a significant 
fraction of the optically dim population at luminosities Lboi 5, 
10 45 e r gs" 1 jMoran et alJl2002l:lHao et alJl2005bl;ISevergnini etal 



veys find that AGN preferentially lie in massive early type hosts 
(references above). As it is well-established that galaxies of dif- 
ferent masses and morphologies cluster differently, these same ef- 



2003 1: iGeorgantopoulos & Georgakakisl I2005f: iGeorgakakis et al . 
2006l : lEckart et alj2006l : lGarcet et alj2007l : lGruppioni et alj2008fh 
Specifically, a large fraction of the sources identified as Type 2 
on the basis of a lack of observed broad lines or optical/IR AGN 
continuum do not show si milar evidence for obscuration based 
on X - ray followup ( see e.g. lAkvlas et al J 120061 : ICaccianiga et al.1 
2008; Hasinger 2008). In many optical/IR studies the Type 2 frac- 
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tion remains as high as /obscured ~ 0.6 — 0.7 up to luminosities 
Z*ol < 3 x 10 45 ergs _1 (L2-iokcV < 10 44 ergs~'), by which point 
the obscured fraction s uggested by X-ray followup has dropped to 

al j |2007n ; iDella Ceca eTait 020041) 



■ 0.2 (see also 



Ta 



er et ; 



/obsci 

and lCaccianiga et al] d2007 ), for example, find that this corresponds 
closely to what is seen in deeper follow-up of X-ray selected 
and apparently op tically dim sources in the XMM-Newton survey. 
lHao et al, (2005b) estimate from a narrow-line (SDSS galaxy) sam- 
ple that the true obscured fraction remains a constant ~ 30 — 40% 
from luminosities L ho \ ~ 10 42 - 10 45 erg s" 1 (Lorn ~ 10 5 - 10 8 L ). 
Thus ~ 50 — 80% of the claimed obscured sources in this lumi- 
nosity range (above Lboi ~ 10 44 ergs _1 orL2_iokeV ~ 10 42 ergs _1 ) 
appear to simply be diluted. 

Note that, although we have found that the obscured fraction 
could be as low as ~ 20%, we are not strongly advocating this 
specific number; this low constant level of obscuration (relatively 
independent of luminosity and redshift) corresponds to the maxi- 
mum effect of dilution on observational surveys, given the uncer- 
tainties in our approach and the observations. The observational 
studies above suggest that the true baseline obscured fraction may 
be ~ 30 — 40% or even slightly higher. Our key point is that the 
observed trend of obscuration with luminosity is very sensitive to 
dilution. And although it has been well appreciated that dilution can 
significantly affect low luminosity AGN samples selected in the op- 
tical/IR, it has not been as widely appreciated that it also biases the 
inferred Eddington ratios, host mass and morphology distributions, 
and even clustering amplitudes (as we have demonstrated). More- 
over it is not clear whether any trivial "correction" can be applied 
to these samples to correct for these effects. 

More detailed observations probing the type of systems we 
appeal to here are clearly important. Unfortunately, they are also 
difficult. As noted above, the sources of interest - those that 
may dominate the "obscured'Vhard population at low luminosi- 
ties in large X-ray surveys - are not analogous to well-studied 
local Seyfert 2 galaxies (which are indeed obscured by dust and 
gas). The latter are primarily nearby, relatively active systems in 
low-mass, star-forming and/or disky host galaxies. In contrast, the 
optically dim/hard population that dominates at low luminosities 
in large X-ray surveys contains a large fraction of massive ellip- 
tical galaxies (M, > 1O H M , and BH mass ~ 10 s - 1O 9 M ) 
firmly in the "red sequence" (i.e. without star formation or signif- 
icant gas content), with low Eddington ratios A ~ 10 -3 — 10~ 2 , 
and relatively high c l ustering amplitudes (se e e.g. Figures 121 121 



and iFiore et al J 120031: iKauffmann et alj 20031; Brand et al 



iNandra et alj 120071: ISilverman et alj 12008 j ; Iffickox et alj 120090 . 



In some ways this population is more analogous to radio galax- 
ies (but many of the systems of interest are still radio-quiet), for 
which a number of models and observations suggest that low Ed- 
dington ratios and differen ces in accretion physics (relative to lo- 



cal S e yferts) are importa nt l Smith et al]|l986r, Falcke & 
19961: lMeieJ l200ll; IhoI [20021; iMaccarone et alj 12001; 



■ & Bicrmann 



Best ct al 



2007; Fender et al.ll2007l : lMerloni & Heinzll2008 

A Lboi ~ 10 44 ergs~' BH in the radiatively inefficient regime 
(with A = Lboi/^Edd ~ 0.001 for m < 0.01) implies a > 1O 9 M 
BH, and there are only a few such objects within ~ 100 Mpc. In- 
terestingly, of the six such objects in the local volume with well- 
meas u red BH masses (fromlTremaine e t al. 2002; Mar coni & Hunt! 
|2003| ; lHaring & Rixl |2004|) . at least four (M87, NGC 1399, 
NGC 4649, and IC 1459) are considered RIAF candidates (al- 
beit with some debate; and in e.g. the c ase of M87 the pres- 
ence of a jet c ompli c ates comparisons; [Sambruna et al. 1999; 
|Pi Matteo etaT] 120001 : lAllen et alj 120001 : ISulkanen & BregmarJ 



l200ll ; iFabbiano et al.|[2003h : they have luminosities L x ~ 10 40 - 
10 43 ergs _1 (Lboi ~ 10 41 - 3 x lO^ergs -1 ) and extremely hard, 
power law-like X-ray SEDs (effective photon indices F — 0.7 — 
1.4). These systems are distinct from many better-studied nearby 
Seyferts in low-mass galaxies, which have softer X-ray spectra V ~ 
2 and Eddington ratios ~ 10~ 2 . These are small number statistics, 
but we reach similar conclusions if we consider the ~ dozen nearby 
BHs down to ~few x 10 s M (se e lEvans et alj|2005l: ISoria et alj 
l2006allbt Imiis & O'Sullivar]|2006l : iGonzalez-Martfn et alj |2006|) . 
Generalizing from this very small sample, if ~ 1/2 of Mbh > 
1O 5 M BHs are active in this luminosity range with X-ray SEDs as 
hard as observed in this subset, this is already sufficient to account 
for > 50% of the observed hard X-ray (apparently "obscured") pop- 
ulation at these luminosities in large samples from z ~ — 1. But 
deep studies searching for e.g. polarized thin-disk signatures could 
set clear lower limits on the fraction of such sources with "tradi- 
tional" obscuration. Clearly, more detailed observations of the (ad- 
mittedly small number of) local analogues is critical to inform our 
speculation that these systems may contribute significantly in large 
surveys. 

There has been considerable debate about the exact physics 
dominant at low Eddington ratios (e.g., an ADAF, a jet- 
dominated spectrum, or some other RIAF model; see, e.g., 



Loewenstein et al. 
Gliozzi et alj|2007l; 



Guainazzi et al. 2000; Angelini et al. 
200 ll; iLumsden et alj|2004 iGallo et aljl2006l; 
Ghosh et alj |2007l : iBrightman & Nandr dl2008h . We stress that we 
are not attempting to advocate any specific model among these. 
Rather, the key question is to test is whether the lack of optical 
emission and the hard X-ray spectra are intrinsic or the result of a 
large column, or even whether or not the X-ray data supports the 
optical conclusions (or whether those can be entirely accounted for 
by dilution). To do so, more observations to break the degenera- 
cies between accretion state, dilution, and genuine obscuration are 
clearly important. Ideally, high-resolution nuclear optical and X- 
ray spectroscopy and optical polarimetry could be applied, as in 
the most well-studied local Seyferts, to break the degeneracies be- 
tween these possibilities and the most popular alternatives ("reced- 
ing torus" models, in which a circumnuclear torus leads to geomet- 
ric obscuration at all luminosities, but with a luminosity-dependent 
covering angle). Although these observational tests may be possi- 
ble for the local analogues noted above, most systems at these low 
luminosities in large surveys cannot be studied in such detail. Ob- 
servations must instead rely on indirect probes; these probes do, 
however, allow statistical tests that are not possible in the small lo- 
cal populations, and the app lication of integra l arguments following 
BH demographics (see e.g. lBrand et al . 2005). 

Perhaps the largest theoretical uncertainty in estimating the 
importance of radiatively inefficient accretion for large AGN sam- 
ples is the rapidity of the transition from a thin disk to a RIAF. If 
the luminosity declines too rapidly in all bands, say, then all ob- 
jects below some rh will simply drop out of all samples, and will 
not introduce significant differences between AGN samples in var- 
ious wavelengths. If, on the other hand, the thin/thick transition 
occurs very slowly with decreasing m, then the SEDs will remain 
relatively X-ray soft and optically/IR bright at lower luminosities 
(Fig. |7J, leading to less of a contribution from such sources to the 
X-ray hard/optically dim populations. Constraints on the contribu- 
tion of such sources to faint "obscured" AGN samples, therefore, 
are a potentially new and powerful way to determine in detail how 
(and indeed whether or not) AGN behavior changes at low accre- 
tion rates. 

If dilution and radiative inefficiency are indeed significant 
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at the levels estimated here, they can naturally account for ob- 
served differences in AGN populations selected at various wave- 
lengths and via different techniques. The apparently contradictory 
host galaxy properties, clustering, and Eddington ratio distributions 
from different samples can, in fact, be reconciled and understood 
using a coherent theoretical framework. The importance these ef- 
fects for large AGN samples is also intimately related to the prob- 
lem of AGN fueling. Black holes that gained most of their mass 
in early, high-accretion rate episodes and have since decayed to 
low luminosities over a Hubble time - typical in merger-driven fu- 
eling scenarios - may show up as low-m populations in massive, 
early-type hosts. On the other hand, systems that continuously or 
stochastically fueled - expected in scenarios where BHs are fueled 
at lower levels by e.g. stellar winds or bars in disk galaxies - will 
plausibly show up as higher-m populations in intermediate mass, 
later-type hosts. 

The results highlighted here also represent an important re- 
vision and caveat to the unified model of AGN. Even if AGN 
physics are perfectly self-similar, dilution implies that objects se- 
lected in optical, infrared, or X-ray samples are not random sub- 
sets of the same parent population (with just viewing angle differ- 
ences). Rather, dilution leads to host galaxy-dependent selection 
effects, biasing e.g. the Eddington ratio distribution, host luminos- 
ity/mass/age/clustering distributions, and other properties. More- 
over, radiatively inefficient populations could represent a large frac- 
tion of the "obscured" population at low luminosities. This should 
not be viewed as a limitation of these observations - rather, it 
presents the opportunity for multiwavelength comparisons to de- 
velop new constraints on AGN fueling mechanisms, feedback, and 
host populations. Since AGN selection may not be strictly inde- 
pendent of host properties at low luminosities, comparison of the 
demographics selected in different wavelengths can break degen- 
eracies between different fueling and host evolution models. 
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